Three Mach Zehnder Interferometers Setup for Production and Observation of 

GHZ-entanglement of Electrons 
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We propose a new single-step scheme for the generation of a GHZ entangled state of three single- 
electron excitations (flying qubits). We also present a method to get a generalized GHZ-state. Our 
idea relies upon the most recent progress in the field of on-demand electron sources and mesoscopic 
Mach-Zehnder interferometry. We also provide the recipe for the unambiguous detection of this 
GHZ state via correlations measurements at the output, which imply the violation of a Bell-type 
inequality which is generalized to the case of three particles. We explain how such measurements 
can be achieved in the context of Mach-Zehnder interferometry, and draw an actual prototype device 
which could be achieved with point contacts and metallic gates placed on a GaAs sample, in the 
integer quantum Hall effect regime. 

PACS numbers: 03.65.Ud, 73.23.-b 



I. INTRODUCTION 

Quantum entanglement, first noted by Einstein- 
Podolsky-Rosen(EPR)i and Schrodinger— is a genuine 
property of quantum mechanics. Qualitative embodi- 
ment of this property was given by Bell 3 - who showed that 
entangled states have stronger correlations than allowed 
by local hidden variable theories (LHVT). Later Clauser 
and coworkers 4 - suggested a more transparent inequal- 
ity (B-CHSH), the violation of which was experimentally 
demonstrated for photons 

During the last decade, several works have been 
achieved in the context of mesoscopic physics to explore 
two particle entanglement. The initial proposals used su- 
perconducting sources of electrons connected to two nor- 
mal metal leads, through which the two constituent elec- 
trons originating from a Cooper pair were split to form 
a singlet state outside the superconductor^. Other pro- 
posals for two electrons entanglement were subsequently 
made using ballistic electrons and point contacts placed 
in the integer quantum Hall effect (IQHE) regime^. Bell 
inequality tests based on stationary current noise cross 
correlations revealed that a maximal violation could be 
achieved £r— 

In the context of quantum optics, Greenberger, Home 
and Zeilinger introduced a maximally entangled tripar- 
tite state, commonly called the GHZ-stateJ^ The Bell 
parameter in the tripartite case can take values up to 4, 
while the corresponding parameter in the bipartite case 
is no more than 2^2^ For both cases the LHVT limit 
remains 2. In this sense the entanglement is considered 
to be "stronger" than for the two particle state, and it is 
thus potentially easier for experimental detection. Ref. 
[HI pointed out that the measurement results with this 
state are in conflict with local realism when quantum 
mechanics makes definite (non statistical) predictions, in 
contrast with EPR states. The GHZ state has a num- 



ber of potential applications in quantum information. It 
can be used for quantum error correction-^, and it was 
proved that tripartite states have advantages compared 
to bipartite ones in quantum teleportation 1 ^ and in dense 
coding 1 ^. There are two established methods to produce 
entangled states of electrons. The first method exploits 
the indistinguishability of fermions and relies on postse- 
lection to get the desired entangled stated (see also Ref. 
I2G for the discussion of the role of projection in this case). 
The second way is to use the interaction between particles 
(implication of this method with flying electronic qubits 
is described in Refs. |21||22I ). The latter has the advantage 
that the evolution during the preparation of the state is 
unitary and thus can be deterministic with up to 100% 
efficiency for producing the desired state. 

While there exist various ways to create and manipu- 
late entangled states of qubits in optical systems^ ' 23 : 24 
and NMR (nuclear magnetic resonance) experiments, 25 
in mesoscopic systems the generation of an experimen- 
tally determined entangled state of few electrons and 
the subsequent proof that Bell inequalities have success- 
fully violated both represent a considerable challenge. A 
proposal which follows faithfully the quantum optics ex- 
perimental protocol^ uses edge states in the IQHE and 
achieves the GHZ state using postselection. 

The purpose of the present work is to go beyond this 
existing protocol and thus to explore three particle or- 
bital entanglement with ballistic electrons propagating 
in mesoscopic devices using the interaction between elec- 
trons. This proposal is motivated by the recent progress 
in achieving single electron sources 2 ^ and for buiding ef- 
fective electron Mach-Zehnder interferometers (MZIs) 27 . 

The outline of the paper is as follows. In Sec. [Hi we 
describe the 3 MZI setup and specify how the beam split- 
ters (BS) and Coulomb interaction operate. In Sec. IIIII 
we justify that a GHZ state is produced at the outcome. 
In Sec. IIVI we describe the Bell inequality test which is 




FIG. 1: (Color online). Proposed setup for the creation of 
an electronic GHZ-state. 3 MZIs interferometers (full lines), 
coupled by Coulomb interactions (ir boxes). Three electrons 
are injected simultaneously in the three MZIs, and propagate 
from the bottom to the top of the figure. The crossing of 
the lines indicate a beam splitter (BS), where an electron can 
exit in one of the two output branches. The dotted lines 
for each MZI show a third BS, used to explain theoretically 
the detection process. The dashed horizontal lines specify 
different levels of the setup which are used in the theoretical 
explanations. 



used to show that maximal entanglement is achieved. In 
Sec. |V]we implement this test for our 3 MZI device. In 
Sec. IVII we show a possible realization of our setup in the 
IQHE regime. We conclude in Sec. IVIII 



II. DESCRIPTION OF THE SETUP 

A schematic drawing of the device is depicted in Fig. [U 
It consists of three MZIs, which are placed side by side 
(shown with full lines in Fig. [IJ. Each MZI consists of 
two incoming channels, which meet at a first beam split- 
ter (BS). The two outgoing channels from these BS prop- 
agate and meet at the second BS. The theoretical setup 
also has a second loop (shown with dotted lines in Fig.[TJ, 
with propagation from the output channels of the second 
BS, and recombination at a third BS. We will use this 
second loop to explain theoretically the detection of the 
GHZ states, and we will show in Sec.|V]that it is possible 
in practice to get rid of this second loop and perform the 
GHZ state production and detection with a single loop. 

We distinguish the right (R) and left (L) side of each 
interferometer and label single electron wave functions 
accordingly. For instance, at any stage of the wave packet 
evolution, ipjc denotes an electron wave packet on the C 
side (C = R, L) of the j-th MZI (j = 1, 2, 3). 

The most likely candidates for electron channels are 
edge states in the IQHE regime, which have the advan- 
tage that they are immune to backscattering effects by 
impurities and have long phase breaking time. Several 



FIG. 2: Reflectionless beam-splitter with transition ampli- 
tudes. 

experiments involving MZI in the IQHE have already 
been performed^!, some of them involving setups with 
two MZIs. A proposal for a detailed setup is given in 
Sec. I VII We assume that a single electron wavepacket 
is emitted in each MZI above the Fermi sea. Such sin- 
gle electron emission was recently experimentally demon- 
strated with the on demand single electron sourc o 26 i 28 
which uses the mesoscopic capacitor as the injector. The 
BS are assumed to be reflectionless, i.e. incident parti- 
cles cannot be backscattered in the same channel, they 
can only be transmitted further "up" (in Fig. [T]). The 
BS can be parametrized by a transfer matrix 

j, _ ftRR t LR \ _ (i sin 6 cos 6< \ , . 

\tRL t LL J \cos9 ismOJ ^ ' 

which relates incoming states to outgoing states and 9 
is the transparency parameter (see Figj^J). For example 
if 6 = the beam-splitter is transparent meaning that 
the incident particle goes from R to L and from L to R 
without scattering, while for 6 = ir/4 the incident par- 
ticle may appear in each of output channels with equal 
probability 1/2. There exists an additional freedom in 
choosing the phases of scattering matrix elements but it 
does not affect results in a crucial manner (see Appendix 
[B]for details). Here we specify the bottom beam-splitters 
to have equal probability 1/2 for transmission in the R 
and L side channels (9 = 7r/4). The transfer-matrix is 
then simply: 

_ 1/V2\ , , 

The main steps in the production of the GHZ-state are 
as follows. The first stage is the synchronized injection 
of electrons (dots with arrows at the bottom of Fig. [TJ, 
which then pass through the first BS of each of the 3 
MZIs. Beyond the first BS, free propagation occurs for 
electrons wave packets, until they reach the "interaction 
region" : the L channel of the first MZ and the R chan- 
nel of the second MZ, and separately the L channel of 
the second MZ and the R channel of the third MZI are 
put in close proximity so that Coulomb interaction effects 
between the two pairs of neighbouring channels become 
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FIG. 3: Illustration of the analogy between a MZI and a spin- 
counter: one MZI coupled to a lead via Coulomb interaction 
(see text for details). 

important. The effective length of this interaction re- 
gion is such that it generates an overall ir phase shift of 
the two particle wave function associated with the two 
neighbouring channels. After the interaction region the 
three electrons propagate freely towards the second BS, 
but we allow to insert a phase difference $ in each MZ 
loop. The phase differences which are needed here could 
be achieved via the Aharonov-Bohm effect, via a scalar 
potential or other means. This phase effectively changes 
the transfer matrix of the second BS. The final stage of 
our setup is the measurement of the state of the three 
particle in the output arms (top of Fig. [1]). 

To understand the role of the interaction regions, let 
us consider first a given MZI with an electron injected 
into the right lead (see Fig. [3]). The electron in this in- 
terferometer could be treated as a "flying qubit", i.e. a 
two-state system, where states are specified by the side 
chosen in interferometer, where the electron travels. The 
initial state therefore is ipR. After the bottom beam- 
splitter, at the level "lvl" (see Fig. [3]), the state of the 
electron is 

* = -j=(i1>R + rl>L). (3) 

Next, we add a wire coupled to the interferometer with 
an electron injected simultaneously. So the initial state 
is ipiR^Pz- When two electrons pass the interaction zone 
they accumulate mutual phase n. After interaction the 
two-particle state remains separable, and each particle 
has a wave-function. If the interaction occurs we get at 
the level "lv2" for the left electron 

Ipir = —^(-ilplR + (4) 

If we do not inject an electron into the right lead then at 
the level "lv2" we have for the left electron 

V^O = -7=(#1.R + V'lL) • (5) 




FIG. 4: (Color online). The 3 MZI setup (up to lv2) repre- 
sented in terms of spin-counters: because of the interaction 
regions, the left and right MZIs in Fig.[T]act like spin-counters 
(large arrows on the figure) coupled to left and right wire of 
middle MZI. 



These two outcomes are orthogonal to each other 
(■0o|Vv) = thus we can distinguish cases when or 
1 electrons travel in the neighbouring lead. Moreover the 
relation (ipol^pTr) = is universal in the sense that it is 
preserved when we change the lead were electrons are 
injected in the MZI or when we place the wire to the 
left of the MZI (as long as the bottom beam-splitter re- 
mains half- reflecting) . We may put into correspondence 
to the wave function <pQ a qubit state | ff) and to (j)^ a 
state | JJ.) . This makes our three MZI setup very similar 
to the situation displayed in Fig. 2J There we have a 
spin-based electron counter with spin flipping after one 
electron passes through a wire (like it is described in Ref. 
[29h . This analogy between MZI and spin- based counter 
described here is helpful for understanding the nature of 
the GHZ-state in the proposed setup (Fig. [I). 

Up to the level "lv2" on Fig. [TJ the evolution can be 
represented in terms of spin-counters (Fig. 2]) After pass- 
ing through the beam-splitter the electron state in the 
middle MZI is if) = -^=(^p R + Depending on the 

arm from which electron came the corresponding spin 
flips. This produces the final state 

¥=-^(|1Wfl|^)+i|Wi|ft>), ( 6 ) 
which is obviously a GHZ-type state. 



III. GHZ STATE PRODUCTION 
A. Standard GHZ-state production 

We give here the detailed explanation about the GHZ- 
state production. We find an exact expression for the 
electron state at the level "lv2" on Fig.[TJ and we identify 
the local unitary transform making it a standard GHZ 
state. "Local" means that the unitary transform is a 
direct product of three unitary operators each acting over 
a Hilbert space of a corresponding electron U = Ui®U2® 

u 3 . 
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Specifically the initial three particle state is chosen to 
be (see bottom of fig. HJ: 

1p — ^1R^2L^3L- (7) 

Note that we choose a simple product state, instead of 
choosing a Slater determinant of the three single particle 
wave functions. This is due to the fact that electron wave 
function do not have mutual parts of trajectories, so that 
exchange effects do not play any considerable role. 

At the level "lvl" after passing the first row of beam- 
splitters: 

* = 2^72 + ^IflXV^fl + #2i)(^3ii + #3i) (8) 

In order to find the state after the Coulomb interaction 
has acted, we rewrite the previous equation as 

1 



2V2 
1 

2yi 



Oil + #i.r)V , 2-rW'3-R + #3l) + 

(tplL+ilplR)ilp2L(lp3H+i''p3L)- (9) 



The state after interaction (at the lv2): 
1 



2V2 
1 



(iplL ~ ilplR)ilp2L(fp3R + itp3L), (10) 



2V2 

or in a more symmetrical way 

* = — 7={i>\L + il/>lR,){—i)lp2R(ilp3R + lp3l) + 



2y/2 
1 

2V2 



{ilplL + 1plR)tp2L(lp3R + #3L), (H) 



Now it is easy to see the unitary transformation which 
produces the canonical GHZ-state 

^ GHZ = -^{lplRlp2Rlp3R + ^\hi>2Li>3h)-, (12) 

from It should transform A=(ipiL + ^lr) into tpiR, 
■^(ipiL-tyiR) into ipi L , -iip2R into ip 2R , ^2L into ip 2L , 
(4>3R - «V>3i) into ip 3 R, and -^(4>3R + #3i) into ip 3L . 
This transform is 



U = E/i ® E7 2 ® t/ 3 , 



(13) 



where 



U 2 



I 



J_ /-» 1 

V2 V 1 ' 
i 
1 



(14) 



So, at the level "lv2" of our setup (Fig. [T]) we have 
the GHZ-state up to a change of basis achieved by the 
transform Eq. 1141 This transform could be performed by 
the second BS row between level "lv2" and level "(lv2')". 

It is interesting to note that from the quantum compu- 
tation point of view we have an algorithm implemented 
on 3 flying qubits which is essentually similar to the cod- 
ing scheme depicted on Fig. 10.2. in Nielsen-Chuangi&. 



B. Generalized GHZ-state production 

In this work we stay focused on the creation of a GHZ- 
state in its original form = -4=(| ttt) + I ill))- How- 
ever in some cases it might be useful to obtain more gen- 
eral GHZ-type state such as |\ff) = a\ ttt) + b \ U4) wit h 
arbitrary a and b. Such a generalized state could be use- 
ful e.g. for quantum error correction^ In our setup it is 
possible to get such a GHZ-type state by a slight modi- 
fication of the above presented method. We simply need 
to replace the half-reflecting BS in the second MZI with 
another one. As we pointed earlier in the general case the 
transfer-matrix of the BS can be represented as follows 



T 



i sin 9 cos 9 
cos 9 i sin 6 



(15) 



With this beam splitter in the second MZI the electronic 
wavefunction is: 

^ = cos9\ t)ib R \ V) +ism9\ ft)), (16) 

which is a generalized GHZ state in a "rotated" basis. 
Using the same unitary transformation U = Ui®U 2 ®Uz, 
where Ui are defined in Eq. (Tl4|) . produces the general- 
ized GHZ-state in its canonical form. 



IV. BELL-TYPE INEQUALITY FOR THREE 
PARTICLES. 



To detect the GHZ-state experimentally we suggest 
to use a Bell-type inequality violation test for three 
particles^. 

Let us discuss some basic facts about Bell-type inequal- 
ities for 3 particles which constitute one of the possible 
generalizations of Bell inequality on a tripartite case. 

We start from an algebraic inequality 

\B\ = 1x^X2X3 + X\x' 2 xz + xiX2x' 3 — x\x' 2 x' 3 \ < 2, (17) 

which is satisfied when Xi,X2,X3, x' ll x' 2 , x' 3 are real vari- 
ables with absolute values less or equal to 1 (this is a 
three particle generalization of the algebraic inequality 
used for Bell inequality tests on two particles). 

Consider a three particle entangled state, written in 
pseudospin notation (R =t, L =4,, each particle can 
be detected in one of two leads L and R). The pro- 
jection of the pseudospin of j on some vector aj cor- 
responds to Xj while the projection on a- corresponds 
to x'j. From the point of view of local hidden variable 
theories, after the creation of the three particle state, be- 
tween the two measurements the first electron has pro- 
jections er ai , cr a ' , the second one a &2 , cr a ^ , and the third 
one CTa 3 , <J & ' 3 (u ai = 1 if the spin is parallel to ai , and 
cr ai = —1 if it has opposite direction). Different out- 
comes of measurements are due to a hidden variable £ 
which varies from measurement to measurement. We 
have identified the real numbers x\, x[, x 2l x 2 ,x 3 , x' 3 with 
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^ai(0> CT ai(£),^a 2 (0' CT ai(£),CTa 3 (£) 5 ^a:,(£)- The aver " 

age value of B is then given by: 

+ cr^cr^cr^ - ^0^0^) p(£)d£, (18) 

where is the distribution function of the hidden vari- 
able. Experimentally this value can be measured by the 
following procedure. Let us define the correlator 

£(a!a 2 a 3 ) = (<r ai a a2 a a3 ) = J (T ai a a2 a aa p{£,)d£- (19) 

Then 

B = E(a' 1 ,a 2 ,a 3 ) + £(a 1; a' 2 , a 3 ) + £(ai, a 2 , a' 3 ) 

-B(ai,a4,4). (20) 

and from Eq. (fTT|) the local hidden variable average is 
such that: 



B < 2. 



(21) 



From a quantum mechanical point of view, the correlator 
of Eq. (|T9"f is an average over the state of an operator 
describing the measurement of three spins in specific di- 
rections: 



and 



^(ai^, a 3 ) = (<7 ai <g><7 a2 ®<7 a 



B = (B), 



(22) 



(23) 



where the analog of the "Bell operator" for three particles 
reads 



B = <7 ai <g) <T a2 ® 0- a3 + (7 ai (g) (T a2 <g (Tag 

+ <T ai <g (T a2 (g - <T ai <g (T a2 (g) a a ' 3 . (24) 

In Appendix [X] we recall under which conditions this 
Bell operator gives a maximal violation of the Bell type 
inequality. The demonstration is based on the fact that 
first, in order for B 2 to have an eigenvalue equal to 16, 
the state must be a linear superposition of | ^tt) an d 
I Hi). Second, the Bell operator B has eigenvalues ±4 
only if the state is related to the GHZ state: 



* = — (lttt> + \m)) 



(25) 



by local unitary transformations. Finally, in this ap- 
pendix, we specify which angles of spin measurement give 
the maximal value for B. The resulting angles are: 



ai 



a 2 



'cost 
sine 







/cos(0 ai ± tt/2) n 
^ = I sin(^ ai ± tt/2) 

'cos(t/> a2 ± tt/2) n 



±tt/2) (26) 







' cos 4> e 
a 3 = I sm(j) a3 



a 3 = I sin(0 a3 ± tt/2) | . 



where ai +4 > a 2 +4>a. 3 = z F 7 t/2 (different signs identify two 
classes of angles corresponding to the upper and lower 
signs in the formulae) . The origin of these classes lies in 
the symmetry of the GHZ-state with respect to reflection 
in the x-y plane. These angles will further be transformed 
into the corresponding beam-splitter parameters 9 and $. 



V. DETECTION SCHEME 
A. Detailed detection scheme 

We are now in a position to describe all the steps for 
production and detection, which are achieved in a rather 
complicated three MZI setup with each MZI containing 
double loops (fig. [J). 

The incoming state of Eq. ([7]) is injected into the first 
BS and becomes that of Eq. © at level "lvl" . The it 
shift is applied to neighbouring channels and it results in 
the rotated GHZ state of Eq. {TTJ at level "lv2" . The 
true GHZ state of Eq. ([12]) is achieved at level "(lv2')" 
by passing through the second beam splitter row (with 
phase shift). The next MZI loop is associated with the 
Bell measurements process. 

To produce the Bell measurement correlator we need 
three values £1,2:2, £3 corresponding to separate mea- 
surements with results within the band [— 1;1]. In the 
spin case, the spin projection measurements serve this 
purpose. In our setup we can detect a particle in the left 
or in the right arm. So one of the possible assignments is 
x = — 1 for the particle detected in left arm and x = 1 for 
particle detected in right arm. Then the Bell correlator 
is an average of a product X1X2X3 which by definition is 



E 



(X!X2X 3 ) =l*Pi + (-l)*P_ 1; 



(27) 



where Pi is the probability for X1X2X3 to be equal to 1 
and P—i to -1 correspondingly. X1X2X3 = 1 in four cases. 
First one is x\ = 1, X2 = 1, £3 = 1, which corresponds 
to the case where the three electrons were observed in 
right arms. We define the probability of this event as 
Prrr- Other three possible cases have corresponding 
probabilities Pllr, Plrl, Pllr- So 



Px=P 1 



RRR. 



P, 



LLR 



P, 



LRL 



P, 



RLL ■ 



(28) 



Analogously 

P-i = Plll + Prrl + Prlr + Plrr- (29) 

Finally, 



E = (2:1X2X3) 



Prrr + Prll + Plrl + Pllr 
—Plll — Plrr — Prlr — Prrl- 



(30) 



This correlator measured on level "lv3" of the setup of 
Fig H] corresponds to the (& z <g> a z <g> <r z ) correlator of the 
pseudospin. 
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Measurement of pseudospin in arbitrary direction. 
As we have seen in Sec. IIII1 in order to make Bell 
measurements we should be able to measure pseudospin 
projections in arbitrary directions, say, n. Nevertheless 
it cannot be achieved directly because we only can 
measure the presence of the particles in the right or the 
left lead ("spin up - spin down"). Instead one can do 
an equivalent measurement which gives the same expec- 
tation values. One should transform the state in such 
a way that the direction n of the pseudospin converts 
to z and then perform the standard measure procedure 
Eq. (j3"0J . This is done by upper the beam-splitters, 
and we obtain: (a ni ® a n2 <S) <J n3 ) for the state on level 
"(lv2')". The relation between rotation parameters and 
beam-splitter characteristics is described in Appendix 
|B1 For each Bell correlator the parameters of this loop 
are different. As soon as we have defined the correlators 
^(a' 1 ,a 2 ,a 3 ),i;(ai,a^,a3),i;(ai,a2,a; j ),£ , (ai,a^,a^) 
we are ready to calculate the value of the Bell observable 

B = E(&' 1 ,a 2 ,a 3 ) + E(Ei 1 ,&' 2 ,a 3 )+E(a 1 ,&2,a' 3 ) 

-£7(ai,a4,a4). (31) 

which is the result of the experiment. 

Granted, the fact that we have to include an addi- 
tional loop in each MZI is quite cumbersome and renders 
the implementation of our setup with mesoscopic devices 
more difficult. In the next section we will show that it is 
possible to perform the same task without the additional 
loop. 



B. Condensed detection scheme 

Here, we insist that the use of an additional loop in 
each MZ can be circumvented, provided that we exploit 
the composition of two BS in series. We will call the setup 
where the additional loop has been removed: "condensed 
detection scheme" . Consider one of three interferome- 
ters (Fig. [TJ). After level "lv2" it has two beam-splitters 
each performing its own unitary transformation. Let the 
transformation matrices of the BS located between level 
"lv2" and level "(lv2')" and the BS between "(lv2')" and 
"lv3" be U and V correspondingly. Then total transfor- 
mation UV is also unitary. The idea is to perform this 
total transformation with a single BS, and thus to keep 
the setup as simple as possible (i.e. remove the dotted 
line part of Fig. [TJ. 

To be precise, note that a single beam-splitter cannot 
reproduce fully all unitary transformations, and there is 
the possibility that the transform UV cannot be repro- 
duced in full generality with a single beam-splitter. But 
it can be shown that, from the point of view of mea- 
surement, it is always possible with a single BS to per- 
form a transformation which gives the same measure as 
the total transformation UV. We leave the details of the 
proof of this statement to the appendices (in particular 
appendix [DJ . 



It is interesting to note that in the condensed detec- 
tion setup, there exists a measurement scheme where the 
desired violation B = 4 is achieved using only changes 
of the phase difference $i of the MZI loops, while the 
transparencies of the BS remain constant: 



' 0i 


= 7T/8, 


o[ = 


vr/8, 




= -tt/2, 




tt/2, 


9 2 


= tt/4, 


e' 2 = 


tt/4, 




= 71", 


<H = 


tt/2, 


9 3 


= 7T/8, 


6> = 


tt/8, 


$ 3 


= -tt/2, 




tt/2, 



(32) 



Here we stress that the $ phases presented in the last for- 
mula are "effective", meaning that they depend on the 
practical implementation of the beam-splitters (values in 
the formula are for T given by Eq. (fTJ). In Appendix [Bl 
we show how to adjust the phases in the case of real ex- 
periments. This constitutes the justification for removing 
the additional loop in each MZ. The whole set of measure- 
ment schemes with B = 4 is presented in Appendix [D] 



VI. SETUP IN THE IQHE 

While the setups which we have presented and dis- 
cussed in this work were schematic, we think that a prac- 
tical implementation is in reach with current experimen- 
tal techniques. Using edge states in the integer quantum 
Hall effect regime (IQHE), one can obtain chiral channels 
where backscattering is impossible. Beam splitters can 
be obtained using quantum point contact (QPC) . Several 
experiments realizing one or two electronic MZI have al- 
ready been achieved j^I. 

The actual geometry for a working device, reproduc- 
ing the setup of Fig. [TJ is shown in Fig. [SJ This pic- 
ture represents a 2D electron gas, with electronic gates 
(in gray) along which edge states propagates. The edge 
states meet at quantum at QPC, whose transparency can 
be controlled by gate voltages. Each MZI has two elec- 
tronic sources (noted Si and S|, i = 1,2,3, on Fig. [SJ, 
and two drains collecting the electronic current (Di and 
D'j). Note that by changing gate voltages, it is possi- 
ble to modulate the length of the different paths, and 
thus to control the Aharonov-Bohm phase acquired by 
the electrons during propagation. The two circles on the 
pictures show interaction regions, where two branches of 
two different MZIs are put close to each other, in such a 
way that a n phase shift is produce when electrons are 
present in the two branches. 

In order to have simultaneous electrons injected in the 
three MZIs, the sources should be well controlled single 
electron sources. Such a source has also already been 
experimentally realized in the IQHE j 26 ' 28 
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FIG. 5: (Color online). Geometry of the setup realized with 
edge states of the IQHE. Sj,S'j (j = 1,2,3) are electron 
sources and Dj , Dj are drains. The edge states meet at quan- 
tum point contact, whose transparency can be controlled with 
gate voltages. The circles show the interaction regions be- 
tween neighbouring MZIs. 



VII. CONCLUSION 

In this article we have demonstrated, how to produce 
a generalized GHZ states for 3 flying qubits, formed as 
orbital electron states in the MZI-type geometries. 

Our setup relies on two state of the art devices of 
nanophysics: a) single electron sources which have been 
demonstrated^ and characterized^ experimentally, as 
well as theoretically 3 -!; b) Mach Zehnder interferome- 
ters, which have equally been the object of thorough 
experimental^ and theoretica l 32 ' 33 investigations. In our 
proposal, both devices should be integrated together, 
and moreover three MZI placed in parallel are needed 
to achieve the GHZ state. While it is plausible to think 
that it will be challenging to build a prototype in the 
near future, we judge that it is useful for the mesoscopic 
physics community to be aware that advanced quantum 
information protocols - here the production of a GHZ 
state - can be achieved with electrons. 

At the beginning of the paper we justified that our 
ideal devices generates the GHZ state. Subsequently, we 
provided a detailed explanation of the type of Bell in- 
equalities which need to be used to prove unambiguously 
that we have generated the proper state. The problem 
which we encountered with the actual MZI setup which 



is needed to implement the Bell test is that it requires 
3 more MZI loops than "condensed scheme", thus mak- 
ing the integration of the device even more challenging. 
Fortunately, we provided a "condensed detector scheme" 
were this complication can be circumvented, arguing that 
the operations achieved by two BS in series can be com- 
bined with a single BS. This allowed to actually draw a 
"realistic" device inspired from the so called "air bridge" 
technique which is used in the experiments of Ref. 1271 . 

This work thus belongs to the ongoing effort called 
"electron quantum optics" where experiments and 
paradigms of quantum optics are reproduced with meso- 
scopic physics devices. The great advantage of the 
present way for producing GHZ states relies in the fact 
that no postselection procedure is needed here. The steps 
described in this paper rely on a unitary evolution of the 
initial state, and moreover, we exploit the Coulomb in- 
teraction between electron wave packets in order to gen- 
erate the desired phase shift on the electron wave func- 
tion. This particular feature departs strongly from the 
photon protocols, as photons interact weakly when trav- 
elling in vacuum. We thus state that the use of electron- 
electron interactions in single electron devices schemes 
may open up new possibilities for quantum information 
schemes which are not envisionable for photons. 

Granted the present work has addressed the case of an 
ideal device which is free from dephasing effects. Such 
phenomena are likely to affect the operation of our de- 
vice as they are already known to be present in single 
MZI electronic setups^ Nevertheless our top priority 
is to clarify how to implement an abstract quantum in- 
formation schemes in an actual mesoscopic device. Our 
analysis could be refined by taking dephasing effects into 
account. The authors plan to consider this in the future. 
In case of similar scheme with two MZI it is done in Ref. 
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Appendix A: GHZ maximal violation parameters 

Here we will prove that in the tripartite case B = 4 can 
be achieved only for the GHZ-state, or the states that can 
be reduced to GHZ with some local unitary transforma- 
tions. Next we will find all possible series of measurement 
angles which give B — 4 for the GHZ-state. We use the 
spin-based approach in this section for convenience. 

First we compute the square the Bell operator (|2~41 and 
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find conditions necessary for (B 2 ) = 16. 

B 2 = 4- [o- ai ,0- a /J ® [o-a^CT^] ® 1 

- 1 ® [(7a 2 , <7ai] ® [<7a 3 7 <7 a <J 

- [ffai^ai] ® 1 ® [^a 3 ,0- a ^]. (Al) 

Pauli matrices (z, j, k = x,y, z) obey angular momentum 
commutation relations: 

[ C 7 i ,a J -] = 2ie^' fe C r fe , (A2) 

from which we obtain the commutator for a spin pointing 
along arbitrary vectors c and d 

[c c ,CTd] = 2i|e|cr e , (A3) 

where e = c x d. For our task we will need vectors 
fi = ai x a' l5 f2 = a2 x a 2 , f3 = a3 x a 3 . They allow to 
represent B 2 in the following form: 

B 2 = 4(1 + |fi||f 2 |cr fl ® C7f 2 ® 1 + |f 2 ||f 3 |l® o- f2 ® cr f3 
+ |fi||fskfi®l®o-fc). (A4) 

This operator has a maximum eigenvalue B 2 — 4(1 + 
|fi||f 2 | + |fa||f 8 | + Ifillfsl) < 4(1 + 1 + 1 + 1) = 16. 

For the maximum value to be achieved we require that 
|fi| = |fa| — | f*3 1 = 1, which implies a x _L a' l7 a 2 _L 
a 2 , a3 _L a 3 . The eigenstates VP of this value obey the 
relations 

cr fl ® cr f2 ® 1* = 

1 ® crf 2 ® cTf 3 * = \P, (A5) 
(T fl ® 1 ® crf 3 * = 

It is convenient to choose a basis in each particle's space 
in relation with the corresponding /-vector. For example, 
| f) for the first spin means that | t) = I t) • ln this 
basis the only states which give B 2 = 16 are | ttt) ancl 
| |||) , so in full generality we choose a superposition: 

* = a| ttt> + b\ III) ■ (A6) 

We can now find which state of the form of Eq. (|A6|) 
is an eigenstate of the operator B with eigenvalue B = 4. 
Choosing the basis we associate the z axis of the first 
spin space with vector fi , and similarly for the other two 
spin states. This implies that ai and a! x are orthogonal 
vectors in the x — y plane, and similarly for a 2 , a 2 and 
for a 3 , a^. For specificity we choose ai = x and a' x = y. 
Analogously for the other spin components. Note that 

<Tai<Ta 2 <Tasl ttt) = ^ 0"a 2 ^a 3 I ttt> = ^ I ttt) 

= -CTa'^ai^l ttt) = *| III), 

0"ai fas 0, a 3 I III) = °ai & a' 2 CT a3 I III) = ^ Ca 2 & a' 3 | III) 

= -fTai^CT^I HI) = -t| ttt), (A7) 

thus implying that 

B(a\ ttt) + 6| III)) = -4»6| ttt) + w| III)). (A8) 



From this we easily confirm that when b = ia we have 
the desired eigenstate with B = 4 

# = -^(|ttt) + *1l44>). (A9) 

Also, note that the orthogonal state \P' = -^(| ttt 

) — i\ III)) has an eigenvalue B — — 4. Both states are 
equivalent to the GHZ state up to local unitary trans- 
formations. In conclusion we have shown that each state 
which can violate the Bell-type inequality maximally is 
a GHZ-state up to trivial transformations. 

Now we proceed with the reverse problem. Starting 
from the GHZ-state, let us find all sets of vectors ai, a 2 , 
a3, a'j, a 2 , a 3 , which yield B = 4. We consider a spin 
correlator (<r ni ® cr n2 ® c"n 3 ) for arbitrary unitary vec- 
tors ni, n 2 , 113 with corresponding spherical coordinates 

^ni , 0ni , $n 2 , 0n 2 , $n 3 , ^n 3 • 

The operator for the spin projection on direction n is 
then 



<r n = n x a x + riyCiy + n z a z = 
This allows us to find 



/ cos(9 n sin0 n e -i M 
Vsin6» n e^" -cos(9 n J ' 
(A10) 



<5- ni ® <T„ 2 ® <T n3 I ttt) = 

cos 6 ni cos 6» n2 cos 6> n3 1 ttt) + 
sin0 ni sin0 n2 sin6» n3 e l ^"i + ^+^ 3 )| |||) + ... 

(AH) 

O-ni ® 0"n 2 ® <7n 3 | III) = 

sin(9 ni sin(9 n2 smO^e-^+^'+^l ftf)- 
cos9 ni cos0„ 2 cos0 n3 | HI) + ... 

(A12) 

and the average of the spin correlator in the GHZ state 
is then: 

(o- ni ®6-„ 2 ®o-„ 3 ) = sin0 ni sin6>„ 2 sin0„ 3 cos(^ ni +^ ni +0„ 3 ) 

(A13) 

The simple formula which we have obtained is a justifi- 
cation of the spin-based approach to Bell type inequal- 
ities. The tripartite case is simpler than the bipartite 
state when the question comes to finding the specific an- 
gles giving Tsirel'son bound, since here all the correlators 
should have unit absolute value. This implies 

sin 8 ni = sin 9 n2 = sin 0„ 3 = 1 , (A14) 

and thus we should consider only vectors of the x—y plane 
as candidates for B — 4. The polar angles of these vectors 
ai, a 2 , as^, a 2 , a 3 are defined by the conditions: 

(<7 ai ® CT a2 ® CT a3 ) = 1, 

(<3- ai ® ^a 2 ® <7 a3 ) = 1, 

(<5- ai ® (T a2 ® <Ta 3 ) = 1, 

(<7 ai ® CT &2 ® <Ta 3 ) = ~1, (A15) 
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which in terms of polar angles reads: 



cos^ 
cos( 
cos( 
cos( 



= 1, 

= 1, 

*a 2 + 0a!,) = 1) 

K>+<t> & >) = -1. 



This system of equations has the solution: 

7T 

0a, + 0h, + 0a „ = ± — • 



(A16) 



(A17) 



with cj) a i — a . =F 7r/2, where i = 1,2,3. All the vectors 
lye in the x — y plane. 

This concludes the definition of the angles giving B = 4 
for the GHZ-state. 



Appendix B: Beam-splitter parametrization 

One of essential parts of the device is the BS. The BS is 
a 4-arm scatterer with a special type of scattering matrix 
which allows to split the beam into two parts without 
any reflection. Generally the scattering matrix of 4-arm 
splitter looks like 



S 



(r\\ t 2 i t 31 tn\ 

t\2 f22 t 3 2 ^42 

£l3 £23 ?"33 ^43 

\tu t 24 t 34 ml 



(Bl) 



In this matrix the tij parameters refer to the transmission 
amplitudes from arm j to arm i, and r.n is the reflection 
amplitude in arm i. This matrix obeys only a unitar- 
ity condition. If time-reversal invariance is added a new 



condition appears tij 
matrix is 



tji. The beam-splitter scattering 



/0 Otsi hi' 

t 32 t 42 

hi t 32 

V*41 *42 



(B2) 



For these parameters the following system of equations 
is valid 



T31 


+ T 41 


= 1, 




T31 


+ T 32 


= 1, 




T 32 


+ T 42 


= 1, 




*31*32 + 


*41*42 


= 0. 


(B3) 



Where TV 



\tij\ 2 . For such an equation set it is rather 



easy to find a simple parametrization. Let us assume that 
* 3 i = cos6»e^ 31 ,t 4 i = sin 6>e^ 41 , t 32 = sin 9e i < t ' 32 , t 42 = 
cos^e 1 ^ 42 . Then we substitute all the t- values into last 
equation and get: 



cos v sm Ue 



1031-1032 



sm f cos tie 



1041-2042 



0. (B4) 



For this equation to be valid must be 

Li + </>42 = 7r + 27m, n g Z. (B5) 



"32 



In the present work we focus on symmetrical beam- 
splitters for simplicity. For a symmetrical beam-splitter 



£31 — *42; 

so the phases are 03i = 
simplifies to: 



"42, ^32 



t 32 = t 41 . (B6) 

041 and Eq. (fB5|) 



031 - 032 = tt/2 + 7m, n E Z. 

One of the possible solutions is 

031 = 0, *31 = £42 = cos 6*, 
032 = 7r/2, t 3 2 = *4i = i sin 9. 



(B7) 



(B8) 



This parametrization will be used in further calcula- 
tions. 

We are interested in the properties of electron trans- 
port from bottom to top, so we will need only the values 
of *3i , t 32 , t 44 , t 42 . Now we can introduce the transfer ma- 
trix: 



T 



( tRR thR 



i sm ( 
cos 6 



cost 
i sin ( 



(B9) 



Note that we have changed indexes because of the reflec- 
tionless nature of the BS: tRR — *32,*rz, = hii^hR = 
*3i ; thL = *3i' In order to find the state of particles after 
the BS one must perform the following substitution with 
wave-functions before the beam-splitter: 



ipR — > cos OipL + i sin 9ipR , 
ipL — > cos 9ipR -j- i sin 9ijj l . 



(BIO) 



Since we have chosen a specific parametrisation of the 
BS, we must confirm that if the BS used in the exper- 
iment has a transfer matrix which is different from the 
one in Eq. (|B9[) it will not spoil the whole experiment. 

We consider again the phase relation in the general 
case: 



031 — 032 — 041 + 042 = 7T + 27m, U g Z 

Each phase can be presented as: 

031 = 031 + 27m = 03i, 

032 = 7r/2 + 032, 
041 = 7T/2 + 041, 



(Bll) 



(B12) 



where the phases 03i, 032, 04i, 042 represent the devia- 
tion from the chosen parametrisation. These obey the 
relation: 



"32 



0. 



(B13) 
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a form similar to the spin rotation operator 



A 








0=0 






FIG. 6: For $ = electron goes to left arm. 



All solutions of this equation are taken into account by 
the following representation: 



»31 
^32 

hi 

i>42 



»3, 

^3 + <f>2, 

h, 

f>i + fa ■ 



(B14) 



where fa, 03, 4>a are some real values. 

This answer has an explicit physical meaning. All 
beam-splitters have the same transfer matrix up to the 
phase accumulation in each channel. This phase accumu- 
lation can affect the phase difference in the left and the 
right paths of the MZI. In order to restore the original 
interference pattern we need to adjust the applied phase 
difference. During the experiment it means choosing a 
reference point with $ = for given MZI. This could be 
achieved using the fact that at $ = an electron injected 
in the right arm of the MZI has zero probability to be 
detected in the right arm (fig. G). 

Beam-splitter with phase. 

In our MZI the electrons have two possible paths to 
travel. If there is a difference in the accumulated phase 
between two paths, it will affect the interference between 
them in a subsequent beam-splitter. It is convenient to 
account for this phase difference (which can be generated 
e.g. by applying a magnetic flux through the MZI loop) 
in the transfer matrix of the BS. If the phase lengths 
of the right and left the paths are fa\ and fai, then the 
transfer matrix of the beam-splitter takes the form: 



T 



i sine 
cos 6 



,i<t>R 



cos tie y 
i sm9e l 



(B15) 



Actually the final results depend only on the phase differ- 
ence of the two paths <£> = </>£ — (f>R so the transfer matrix 
can be multiplied by any phase factor without affecting 
the final results. We can represent the transfer matrix in 



T 



sintfe-'*/ 2 -icosfle 4 */ 2 
^cosfle- 1 */ 2 sin^e**/ 2 



(ir-29) 



(B16) 

Since the measurement procedure in the condensed 
scheme includes changing 9, $ and measuring (a z <E)(J Z <Ei 
a z ) at level "(lv2')" (Fig. [T]) it is convenient to define a 
measurement operator with respect to level "lv2" of the 
scheme, where the electron state is invariant during the 
measurement procedure. 

As we see the operator is a sum of two rotations, first 
a rotation of angle $ around the z axis, second another 
rotation of angle n — 29 around the x axis. The rotation 
around x-axis corresponds to the action of the beam- 
splitter itself, while the rotation around z corresponds 
to an Aharonov-Bohm phase accumulation. This whole 
rotation transforms a vector 



'sin 29 sin $ N 
n = | sin 29 cos <E> 

- cos 29 



(B17) 



unto the z axis. As the vector thus defined is in one- 
to-one correspondence with the pair of BS parameters 
9, <J> we call it a characteristic vector of the BS. Rotation 
means that 



Tin) 



t) 



(B18) 



The projection of the rotated pseudospin onto z can be 
shown to be equal to the projection of the initial pseu- 
dospin onto n. 

(U a ^\a z \U n ^) = (¥|E#_>,<t,17»-«|¥) = (tf |<7„|tf). 

(B19) 

For three-particle spin-correlators the same relation 
holds 



{^\T+a z T x ® T+a z T 2 ® T+a z T 3 \V) = 



<7„ 3 |*) (B20) 



So now we have found that the measured correlator 
with respect to level "lv2" is (cr ni (8) cr n2 <g> <J n3 ). 



Appendix C: Pseudospin analogy 

Previously in appendix[X]we have defined angles yield- 
ing B = 4 for three-spin GHZ-state. To obtain the pa- 
rameters of the actual setup (Fig. [T]) we need to "trans- 
late" the spin angles into BS parameters. This chapter 
contains the derivation of the "translation" rules. 

We start with the logic of the GHZ-experiment de- 
scribed in Sec. IIIII We need three values Xi,x 2 ,x 3 corre- 
sponding to separate measurements with results within 
the interval [— 1; 1], In the spin case, the spin projection 
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measurements serve this purpose. In our setup we can 
detect the particle in the left or the right arm. So one of 
possible assignments is x — —1 for the particle detected 
in left arm and x — 1 for particle detected in right arm. 
Then the Bell correlator is 



{xix 2 x z ) = Prrr + Prll + Plrl + Pllr 

—Plll — Plrr — Prlr. — Prrl- (CI) 

This correlator measured on level "lv3" of the 
setup (Fig fT]) corresponds to the (<7 Z <8><7 J! (8>(7 J! ) correlator 
of pseudospin. 

Due to the rotation between "(lv2')" and level "lv3" 
the correlator of Eq. (|C1|) corresponds to (a ni (^a n2 ^a n3 ) 
measured on the initial pseudospin state at "(lv2')" which 
is GHZ. The vectors 111,112,113 are characteristic vectors 
of the corresponding beam-splitters (|B17|) . 

Now we calculate BS matrices corresponding to the 
angles giving maximum violation (|A17|) . All the vectors 
are in the x — y plane, thus 8 — ir/4 for all of them. For 
the vector with the polar angle <f> corresponding $ = a = 
7r/2 — </> and the BS transfer matrix is 




where equivalence means that the overall phase does 
not change the observable outcomes, a parameters then 
should obey the restrictions: 



Q!ai + «a 2 + a a3 = ± 7T/2, 

a a { = a ai ± 7r/2, 



(C3) 



for i = 1,2,3. 



Appendix D: Settings giving maximum violation of 
Bell inequality. 

In appendix [C] we have defined the properties of BS 
between levels "(lv2')" and "lv3" of Fig. □ The next step 
lies in defining the parameters of BS in the condensed 
scheme. For this purpose we use the properties of unitary 
transformation made between levels "lv2" and "(lv2')" of 
the initial setup (fT4|). 

The problem is simplified by the fact that all BS be- 
tween levels "(lv2')" and "lv3" are of the same type 



T(a) 



s/2 \/2 c 

72 71 ( 



(Dl) 



First and third MZI. 

Since U3 = U\ all the results are valid for the third MZI 
as well. The first MZI chain makes two consequent uni- 

f—i 1 

tary transformations, the first one is U\ = I ^ . 



next is T(a) = -h= 



V2 I 1 ie l 



The product of these is 



unitary 

Tx{a) 



T(a)Ui 



6 

- ie % 



1 



(D2) 



Shifting the whole matrix phase by —a/2 we get following 
expression 



cos(a/2) sin(a/2) 
— sin(a/2) cos(a/2) 



(D3) 



The next step is to find from the matrix the beam-splitter 
the parameters 8, $. We encounter the problem that 
the matrix of Eq. (|D3[) cannot represent the transfer- 
matrix of any BS. The relation 8.Yg(tRR/tRL) = ~ which 
is universal for any representation of the transfer matrix 
(see Eq. ([B15]) ) is not valid for Eq. (|D3]) . The BS transfer 
matrix should not necessarily be equal to Ti($) but it 
should lead to the same three-particle probabilities. This 
happens if transfer matrix of beam-splitter equals to: 



Ti(*, 



1) 



cos(a/2)e* sin(a/2)e*\ 
-sin(a/2)e^ 2 cos(a/2)e^ 2 ' ^ ' 



where (f>i and 4>2 are some phase shifts. They do not con- 
tribute to tripartite probabilities while this matrix with 
properly chosen phase shifts can be transfer matrix of 
BS. 

If tan(a/2) > then <\>\ — tt/2, 4>2 — ^ and 



$ = -. 
2 

If tan(a/2) < then fa = tt/2, 4> 2 = and 

7r + a 



(D5) 



$ = . 

2 

Both formulas can be unified as: 

7r — signtan(a/2)a 



$ = signtan(a/2) — . 
If a € [— 7r; 7r] last formula is simplified 



(D6) 



(D7) 



2 

7r|a| 



(D8) 



2q 



Second MZI. 

For second MZI everything is simpler since U2 is a 
diagonal matrix. 



T 2 (a) =T{a)U 2 



V2 



(D9) 



12 



This matrix corresponds to a BS with parameters — 
tt/4 and $ = a - tt/2. 

Finally we can introduce the algorithm for the gener- 
ation of the measurement settings resulting in B = 4. 



1 . Choose angles a ai , a &2 , ot &3 obeying relations 



a ai + a a2 + a &3 = tt/2 ± tt/2, 

"af = ± TT/2, 



(D10) 



where i — 1,2, 3. 



2. With necessary 27r shifts make a a /,a a . be in 



• In the second MZI the transparency of the BS is 
not changing during the experiment; 

Among measurement schemes there exist ones where 
transparencies of all BS do not change during experiment. 
For example, choosing 



(D12) 



we get 



«ai = 


-37T/4, 


a ai =3tt/4, 


«a 2 = 


-tt/2, 


"a!, =71", 


a a3 = 


-37T/4, 


a«i =3tt/4, 



3. Calculate BS parameters 



Pi, 3 — 

$1,3 = 
$ 2 



*ai, 3 | 



2a 
tt/4, 



7 1,3 = ■ 



ai,3 



2 

""l^ai,. 

2av_ . 



tt/2, 



e' 2 = tt/4, 
* 2 = s _ n / 2 > 



(Dll) 



From the last equation we can extract following com- 
mon features of these measurement sets. 

• In the first and the third MZI the AB phase can 
take values ±7r/2; 





= 7T/8, 


*i = 


tt/8, 


$1 


= -tt/2, 


$' 1 = 


tt/2, 




= tt/4, 


02 = 


tt/4, 




= 7T, 


$2 = 


tt/2, 


3 


= 7T/8, 


0.3 = 


tt/8, 


$3 


= -tt/2, 


$3 = 


tt/2, 



(D13) 



In conclusion, we found the general parametrization 
for measurement settings giving B = 4 in a regular MZ 
setup, and described a special case where we get a max- 
imum violation only by adjusting AB phases. 
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